Classification of cosmology with arbitrary matter in the Hofava-Lifshitz model 
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In this work, we discuss the cosmological evolutions in the nonrelativistic and possibly renormal- 
izable gravitational theory, called the Hofava-Lifshitz (HL) theory. We consider the original HL 
model (type I), and the modified version obtained by an analytic continuation of parameters (type 
II). We classify the possible cosmological evolutions with arbitrary matter. We will find a variety of 
cosmology. 
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I. INTRODUCTION 

Studies on ultra-violet (UV) completion of gravity have 
a long history. A lot of understanding has been obtained 
especially from string theory. Very recently, a new class 
of UV complete theory of gravity was proposed by Hof ava 
by generalizing the ideas discussed in Ref. 0. The 
theory does not have the full difFcomorphism invariance 
as in Einstein's general relativity, but only has a local 
Galilean invariance, and Einstein's relativity may emerge 
at the infra-red (IR) fixed point. The model has been 
motivated by a scalar field model discussed by Lifshitz Q , 
to explain the quantum critical phenomena in condensed 
matter physics, in which the action has z = 2, where z 
represents the dynamical critical exponent given by 



t tt, 



(1) 



The model proposed in Ref. [l| has the scaling dimension 
= 3 and is often called the Hofava-Lifshitz (HL) model. 

There would be many implications of the HL model for 
cosmology, which were firstly discussed in d, Q . In Ref. 
, it was suggested that the divergence of speed of light 
in the UV region may resolve the horizon problem. The 
possibility of generation of the scale invariant spectrum 
in the UV regime without inflation was suggested in 0- 
Q . Some detailed studies on perturbations in UV regime 
have been presented in Refs. [H, 0, [ol-flll . Tensor per- 
turbations were investigated in Ref. [l^l, which showed 
that primordial gravitational waves are circularly polar- 
ized due to parity violation. Phenomenological implica- 
tions were suggested in [l^ . The black hole solutions and 
related issues have been studied very actively in [3? - 
\v\ . Some theoretical issues and generalizations of the 
HL model have been discussed in [20l - |24| . 

This article is devoted to discuss the cosmological evo- 
lutions in HL model. We consider the original HL model 
(type I) , and the modified model obtained by an analytic 
continuation of parameters (type H). We classify the evo- 
lutions of the universe with arbitrary kind of matter. The 
article is constructed as follows. In Sec. II, we briefly re- 
view the HL model and give equations of motion in the 
cosmological background. In Sec. HI, we discuss the 
cosmological evolutions in analogy with a point particle 



moving in a potential in classical mechanics. We classify 
the possible cosmological evolutions. In Sec. IV, we give 
a brief summary before closing the article. 



II. HORAVA-LIFSHITZ MODEL AND 
COSMOLOGY 

We briefly review the HL model. The metric which has 
the local Galilean invariance can be written in the form 



-TV^dt^ + {dx^ TV'dij [dx^ ^ NHt^ ,(2) 



and Ni = g^N^ which is similar to the ADM decom- 
position in Einstein's general relativity 25|. The scaling 
dimensions of various quantities in the momentum unit 
are given by [t] = —3, [a;'] = —1, \N\ = \gij\ = and 
[TV'] = 2. The dynamical variables are A^, Ni and 5^-, 
which is very similar to the ADM decomposition in Ein- 
stein's general relativity [2^ . 

The total action is composed of the kinetic part and 
potential one: S = Sk + Sy . The kinetic action is given 

by 



Sk = ^ 



dtd^x 



(3) 



and K''^ = g'^'^g^^K^, where the extrinsic curvature is 
defined by 



1 /d 
m\di 



V.AT,). 



(4) 



V; denotes the covariant derivative with respect to gij. 
Correspondingly, with the detailed-balance condition, 
the potential action is given as 



Sv 



dtd^x^N 



2w4 



'^'^ ''''' R \7 "'^^ 
rUif VjK k 7; — 



2w^ ^ 
8(1 -3A) 



1 -4A 



+AR- 3A2 



(5) 



2 



where Rij is Ricci tensor associated with g^, and the 
Cotton tensor Cij is transverse, conserved and vanishing 
for ah conformally flat spaces and has scaling dimension 
3. The full action is given by 



Si — Sk + Sv 



-z-^—=Rie'^jR^k — -^RijR^^ 
2w^ y/g 8 



8(1 -3A) 



(- 



4A 



-R^ +AR- SA^ 



(6) 



where scaling dimensions of parameters are given by [k] — 
0, [w] = 0, [/z] = 1, and [A] = 2. A represents a dynamical 
coupling constant, susceptible to quantum corrections 
For convenience, one can rewrite the above action 



Si 



where 



dtd^Xy/gN 

ijk 



XK- 



i—RriV,R\ + QRrjR'' + +^R + <y (7) 
v9 



K 

1 -4A 



8(1 - 3A) 



3A) 4 ' 
( - 3A2) 
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2«;2' ^ 



U-3A) 



A, 



(8) 



The action is invariant under the restricted class of dif- 
feomorphisms (foliation-preserving diffeomorphism) t' = 

h{t) and {x'Y = h^t.x"-). For A 1/3, the theory has 
classical, anisotropic conformal invariance. Comparing 
the action with that of general relativity in ADM for- 
malism, one can read the emergent constants 



9 



A 



1-3A' 



G 



327rc' 



A^ = -A, 



(9) 



appearing in IR. Einstein's general relativity is recovered 
for A = 1. 

Following Ref. [3l , one can make an analytic continu- 
ation of parameters /i — > i/x and iv^ — > —i'uP'. Then, the 
above action can be rewritten as 



+ l2^R^iyJRk + C2R^JR'' + mR + 6^ +{m 

where new parameters arc defined by a2 ~ a, P2 = 
-P, 12 = -7, C2 = -C, m = -?7, 6 = and 
(72 = — cr. The emergent speed of light is given by 
c = (K2/i/4)-\/A/ (3A — 1). For convenience, in this ar- 
ticle, we call the model Eq. ([7]) "type I" and the model 
Eq. (Uni) "type IF. 



For the cosmological background, i.e., the homoge- 
neous and isotropic background, N = N{t), Ni = and 
gij = a'^(t)jij, where jij is maximally symmetric 3-space 
whose curvature is given by R]j = 2kjij and R'^ = 6k. In 
the cosmological spacetime, the Cotton tensor vanishes 
and the equations of motion get simplified. 

In the type I theory, the equations of motion (Fricd- 
mann equations) arc given by 



H 



2 KV^(fc-a^A)2 



P 



16(l-3A)2a4 3a(3A-l)' 

{k + 3a'^A){k-a'^A) 



P 



a(3A - 1)' 



(11) 



where H := d/{aN) (in the later discussion, we sim- 
ply set = 1) and p and p are the energy density and 
pressure of the matter sector. In the type II theory, by 
changing 



-/i , the equations of motion become 



H 



2 KV^(fc- 



a^A)^ 



16(l-3A)2a4 3a(3A-l)' 



2[H+-H' 



a(3A- 1) 



\&a\l - 3A)2 



{k + ia'^A){k- a^h) 
(12) 



It is convenient to represent the Fricdmann equation in 
analogy with a dynamics of a point particle in a potential 
in the classical mechanics, -|- Vm{a) ~ Q {m = I, II), 
where 



pa 



AvV(fc-Q2A)2 

16(l-3A)2a2 3a(3A-l) 



, (13) 



where ei = +1 and en = — 1, respectively. In the absence 
of the matter p = 0, the effective potential is not sensitive 
to the value of A unless A = 1/3. With matter, it depends 
on whether A is larger or smaller than 1/3. In the regime 
A > 1/3, since a > 0, the sign of the energy density terms 
in Eq. ([T^ is positive. While, in the regime where A < 
1/3, the sign of the same term is flipped. In the case A = 
1/3, where the theory develops an anisotropic conformal 
invariance, the scale factor becomes non-dynamical as 
seen in Eq. pT|) and P^ . In the later discussions, we 
assume A 7^ 1/3. Note that the case of A < 1/3 induces 
the repulsive gravitational force (see ()13p ). and also the 
perturbation about the flat background provides a ghost- 
like scalar mode [l[, which is potentially dangerous if it 
is coupled to the matter. Thus, this branch may not 
be realistic. But for the mathematical completeness, in 
this article, we will include the case of A < 1/3 into our 
analysis. 



III. COSMOLOGICAL EVOLUTIONS 

In this section, we discuss the classification of the cos- 
mological evolutions with matter in the HL model. We 
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assume that the energy density of the matter is nonncg- 
ative p > 0, which can be parameterized as p = poag/a"-, 
where n is arbitrary integer and po ^ 0. One might think 
that it would be enough to consider the matter of 7i = 6 
with the equation of state p = p, because in the nonrel- 
ativistic theory with the dynamical critical exponent z, 
the equation of state of the massless particles becomes 
p = {z/3)p. But it seems to be too restrictive. The HL 
theory itself is purely a gravitational theory and does not 
specify the nature of the matter. In addition, there is no 
particular reason why the matter sector may not share 
the same dynamical critical exponent with the gravity 
sector, although it may be plausible. Therefore, for com- 
pleteness of our analysis, it is appropriate to keep n to 
be an arbitrary integer. 

It is important to note two important properties. In 
the type I model and in the presence of the matter p > 
there is no solution for A < 1/3, because the potential 
Vi{a) defined by is always positive. In the type II 
model and in the presence of the matter p > there are 
always monotonic solutions ^ , for A > 1/3, because the 
potential Vii{a) defined by is always negative. Thus, 
in the presence of the matter, we will basically focus on 
the case of A > 1/3 for the type I model and on the case 
of A < 1/3 for the type II model. 

In the simplest case fc = A = 0, both in type I and II 
theories, the solutions are given by 



p^n^al 



,2/n 



12a(3A- 1)^ 
which is only possible for A > 1/3. 



(14) 



that the universe starts from a = and ever approaches 
a = k/ A. The second possibility is that the universe 
starts to co ntract from the infinity and ever approaches 
a = yjk/ K. The third and fourth ones are oppositely that 
it starts from a = i/fc/A, and then collapses at a = or 
ever expands toward a — oo, respectively. 

For k = (but A 7^ 0), there is the dc Sitter solution 



aii{t) = ao 



,±Hat 



K^/iA 



4(3A- 1) 



(16) 



For A = (but fc 7^ 0), the universe is just like the 
radiation-dominated one in the relativistic theory, where 
the solution is given by 



aii{t) 



K^pk 



2(3A - 1) 



1/2 



,1/2 



(17) 



The vacuum case 



B. The case n > 5 

Then, we study the cosmology with matter. In this 
section, we discuss the case n > 5. The radiation energy 
density would scale as in the UV regime, not as 
in the relativistic cosmology, because in the UV regime 
the dispersion relation of a massless particle becomes w cx 
p^, where p is the physical momentum and is redshifted 
as p cx 1/a, and the radiation energy density is given by 
p ^ nuj (X 1/a^, where n oc is number density of 

massless particles. 

Type I model (X > 1/3): For all the choices of A and 
fc, big crunch solutions are obtained. 

Type 11 model (X < 1/3 j; For all the choices of A and 
fc, bouncing solutions are obtained. 



Firstly, we consider the case of a vacuum solution p = 

0. 

a. Type I model For the nonzero k and A, there is 
only the static solution for fc/A > 0, which is given by 
ai = y/k/A. where the effective potential Eq. ([13]) is 
vanishing. 

For fc = (but A 7^ 0) or for A = (but k = 0), 
the potential becomes positive everywhere and there is 
no cosmological solution. 

b. Type II model: For k ^ and fc/A > 0, there is 
an exact solution given by 



ajj{t) 



- cie 



±2Hat 



Hn := 



K^pA 



4(3A- 1) 



,(15) 



where ci is an arbitrary constant, which can be positive 
or negative. There are four possibilities. The first one is 



^ By " monotonic" , we mean the universe which expands from a = 
to a = 00 and contracts from a = 00 to a = 0, without any 
collapse or bounce. 



C. The case n = 4 

This case corresponds to radiation contribution in IR 
regime. 

Type I model (X> 1/3): If 



K^p^k"^ ^ poflo 



16(3A-1) 3a 
there is no solution, or there are cyclic universes. If 



< 



Poat 



16(3A-1) 3a ' 



(18) 



(19) 



there are big crunch solutions. 

For fc = (but A 7^ 0), which is the special case of 
Eq. ([T9l) . there are big crunch solutions. For A = (but 
fc 7^ 0), there is no solution for the case Eq. (|18p and 
there are monotonic solutions for the case Eq. ((T9| . 

The exact solution is given by 



aiit) = 



1 



V2H0 



K - (A'2 - 4CiJ2)i/2 cos{2H{t - 



4 



where 



C 



16(1 -3A)2 3a(l-3A) ' 



16(1 -3A)- 



r, A' 



8(1-3A)2 



(21) 



Type II model: For A < 1/3, there are monotonic, 
or big crunch, or bouncing solutions. For fc = (but 
A ^ 0), which just the special case of Eq. (f^. there are 
bouncing solutions. For the case of A = (but fc 7^ 0), 
there arc big crunch solutions. 

For A > 1/3 and = (but A 7^ 0), there is the exact 
solution 



The solution only exists for > 4CH^, which is con- 
sistent with A > 1/3. 

Type II model (X<l/3): If 



aii{t) 



I K^nA V 



K2(3A-l)poag 



X sinh 



16(1 -3A) 3a 
there are bouncing solutions. If 



16(1 -3A) 3a 



(22) 



(23) 



/iA V 6 
|- 3k2^A 
8(3A- 1) 



2/3 



(27) 



there are monotonic solutions. 

For fc = (but A 7^ 0), which is just the special case of 
Eq. ([22]) . there are bouncing solutions. For A = (but 
A: ^ 0), there is no solution for Eq. <{22\ and there are 
monotonic solutions for Eq. ([251) . 

The exact solution is given by 



E. The case n = 1 and n = 2 

Type I model (X > 1/3 j; For A > 1/3, in this case, 
there is no solution, or there are cyclic solutions. For 
fc = (but A 7^ 0), there are big crunch solutions. For 
A = (but fc 7^ 0), there are bouncing solutions. 

Type II model (X < l/3j; There are monotonic, or big 
crunch, or bouncing solutions. For fc = (but A 7^ 0), 
there are bouncing solutions. For A = (but fc 7^ 0), 
there are big crunch solutions. 



au{t) 



1 



for K"^ > AC Hi, and 
1 



K±(K^- ACHiy/^ cosh{2H{t - io)I?4) 



F. The case 7i = 

Type I model (X>l/3): For A 7^ and fc 7^ 0, if 
^V'A2 po 



> 



aii{t) 



V2H0 

for K'^ < AC Hi, where 

K'^p?k' 



16(3A-1) 3a' 

K ± {ACHl - A'2)i/2 smh{2H{t - to)I?5) there is no solution, or there are cychc solutions. If 

h p A Po 



< 



16(3A-1) 3a 



(28) 



(29) 



C 



16(1 -3A)2 3a(3A-l)' 



K^p^A^ 
16(1 -3A): 



n'^p'^kA 
K1-3A)2 



(26) 



there are bouncing solutions. 

For fc = (but A 7^ 0), if Eq. there is no solution, 
and if Eq. ([2^ there are monotonic solutions. For A = 
(but fc 7^ 0), there are bouncing solutions. 

The exact solution is given by Eq. ((2(1)) where 



The first and second solution corresponds to A < 1/3 
(bouncing solutions) and A > 1/3 (monotonic solutions). 



D. The case n = 3 

The case n — i corresponds to that of massive parti- 
cles: p ^ mn cx \/a^, where n is the number density of 
particles with mass m. 

Type I model (X > 1/3/; In the type I model, there 
is no cosmological solution unless A > 1/3. For A < 1/3 
In this case, there is no solution or a cyclic universe is 
available. For fc = (but A 7^ 0), big crunch solutions are 
obtained, while for A = (but fc 7^ 0) bouncing universes 
are. 



c 



H?i:= 



16(1 -3A)2^ 
kV'A2 



K 



n'^p^kA 
8(1-3A)2 

Po 



16(1 -3A)2 3a(l-3A) 



(30) 



The solution only exists for > ACH^, which is con- 
sistent with A > 1/3. 

Type II model (X<l/3): For A 7^ and fc 7^ 0, if 



^V'A2 ^ PO 



16(1 -3A) 3a' 



(31) 



there are monotonic, or big crunch or bouncing solutions. 
If 



«^M^A2 ^ Po 



16(1 -3A) 3a' 



(32) 
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there are big crunch solutions. 

For k = (but A 7^ 0), if Eq. (|3T|) there are monotonic 
solutions, and if Eq. p9|) there is no solution. For A = 
(but k ^ 0), there are big crunch solutions. 

The exact solution is given by Eq. for > 

4:CH§, and Eq. ^ for < iCH^, where 



are big crunch solutions for n > 1. The case n = 
is marginal: If Eq. (|28p there is no solution, but if 
Eq. ([29|l . there are monotonic solutions. 



C 



K 



16(1 - 3A) 
' 8(1 -3A)2 



2' ^0 



16(1 -3A)2 3a(l-3A) ' 
(33) 



The first and second solution corresponds to A < 1/3 and 
A > 1/3. 



G. The case n < — 1 

Type I model (X > 1/3): For all the choices of A and 
fc, there are bouncing solutions. 

Type II model (X < 1/3): For all the choices of A and 
k, there are big crunch solutions 



IV. SUMMARY 

We investigated evolutions of homogeneous and 
isotropic universe in the recently proposed nonrela- 
tivistic, renormalizable gravitational theory with an 
anisotropic scaling z = 3 (called the Hoi^ava-Lifshitz (HL) 
theory), where z is the dynamical critical exponent de- 
fined in Eq. ((U . We considered the original theory (type 
I) and the modified model obtained by an analytic contin- 
uations of parameters in the original theory (type II) . The 
dimensionless, dynamical coupling parameter A is con- 
tained into the theory, which is sensible to the quantum 
corrections. The theory has IR fixed point and for A = 1 
the theory coincides with general relativity in IR regime. 
For A = 1/3, the theory has an anisotropic conformal 
invariance and scale factor becomes non-dynamical. In 
this article, we did not restrict the range of A. We also 
have assumed that the matter energy density is positive 
p>0. 

In the type I model, we found that 

• For the vacuum universe, there is a static solution, 
a = v/fc/A, for A 7^ and fc/A > . 

• With arbitrary kind of matter, there is no cosmo- 
logical solution for A < 1 /3 

• For A = (but k ^ 0) and with matter (A > 1/3), 
there are bouncing solutions for 71 < 3 Q , and there 
are big crunch solutions for 71 > 5. The case n = 4 
is marginal: If Eq. ([T8| there is no solution, but if 
Eq. (fT9|) . there are monotonic solutions. 



• For fc = (but A 7^ 0) and with matter (A > 1/3), 
there arc bouncing solutions for n < — 1, and there 



In the type II model, we foimd that 

• For A 7^ (and fc/A > 0) and in the absence of 
matter, there are exact solutions given by Eq. ([15]) 
[3| • They are solutions approaching a = \/k/ A or 
starting from a = \Jk/ A. 

• With arbitrary kind of matter (A > 1/3), there are 
monotonic solutions for any n. 

• For A = (but fc 7^ 0) and with matter (A < 1/3), 
there are big crunch solutions for < 3 Q, and 
there are bouncing solutions for n > 5. The case 
n = 4 is marginal: If Eq. ([22]) there is no solution, 
but if Eq. ([23| . there are monotonic solutions. 

• For fc = (but A 7^ 0) and with matter (A < 1/3), 
there are big crunch solutions for n < ~1, and there 
are bouncing solutions for n > 1. The case n = 
is marginal: If Eq. ([32]) there is no solution, but if 
Eq. pi|) . there are monotonic solutions. 

It may be helpful to clarify the characteristic differ- 
ences of the cosmological evolutions in the HL theory 
from those in general relativity. 

• In the case of the type I theory: supposing that 
p > and A > 1/3 (the gravity is attractive), in 
the case of fc = and A = 0, there is no essen- 
tial difference from the case of general relativity. 
In the case of nonzero fc, the particular difference 
from the case of general relativity is the existence 
of the bouncing evolutions in the absence of mat- 
ter which breaks the weak energy condition. Note 
that in general relativity, it is impossible to real- 
ize the bouncing Universe under the weak energy 
condition. In the HL theory, the corrections due 
to the higher spatial curvature terms could effec- 
tively break it. Therefore, this is the novel effect 
in the HL gravity. On the other hand, even in the 
flat (fc = 0) or open (fc < 0) Universe, with matter 
of n > 0, the presence of non-zero A induces the 
maximal size of the Universe, i.e., the big crunch 
solutions. 

• In the case of the type II theory: supposing that 
p > and A > 1/3 (the gravity is attractive), in the 
case of fc = and A = 0, there is no essential differ- 
ence from the case of general relativity. However, 
in all the cases, the effective potential Eq. (|13p is 
always non-positive. Thus, even in the case of the 
closed (fc > 0) Universe, there is no recollapsing 
Universe, which is the crucial difference from the 
case of general relativity. 

In this article, for simplicity, the detailed balance con- 
dition was assumed, but it may not be essential. The 
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absence of the detailed balance condition introduces an 
additional contribution to the effective potential Eq. 
(fT3|) (See e.g.,!!^), and leads to a richer variety of the 
cosmological evolutions. In addition, in the type I theory 
with the detailed balance condition, to obtain the attrac- 
tive gravitational force, one has to choose A > 1/3 (see 
Eq. (|13p ). Then in order to have the real emergent speed 
of light, from Eq. ^ the condition of A < has to be 
imposed. Namely there is only the anti-de Sitter vacuum 
{Ae < 0). Although this is not the case for the type II 
theory, the detailed balance may not be suitalbe to obtain 
the realistic cosmology. For the detailed classifications in 
the case without the detailed balance condition, see [1^ . 

Before closing this article, we shall note some proper- 
ties of the model, which would be important for future 
studies. We mainly focused on the possible cosmologi- 
cal evolutions in the UV regime and ignored the process 
of the recovery of the cosmology in Einstein's relativity. 
The UV cosmology should be smoothly matched to that 
of IR. But, there still seem to be some unclear points 
about this connection. First, we assumed that A is a 
constant but actually A would be sensitive to the pertur- 
bative corrections. A would eventually approach some IR 
fixed point and variation of A would affect the cosmology. 
Second, in Ref. P, [3, 0, [111, [iBl , it has been pointed out 



that in the HL theory there is another (scalar) degree of 
freedom other than those appearing in general relativ- 
ity, which corresponds to a scalar mode in the linearized 
theory. The analysis of perturbations about a flat space- 
time indicates that the kinetic term of this scalar mode 
becomes a ghost for A < 1/3 and A > 1, and may be im- 
portant for stability of the solution, if this mode is cou- 
pled to the matter perturbations [J [i3 • It also has 
been pointed out that in the cosmological background 
this mode may be useful to produce the scale invariant 
spectrum without inflaton in the UV regime, if the de- 
tailed balance condition is broken [9|. The decoupling of 
scalar mode in approaching IR would also be essential to 
recover the relativity. 
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